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Abstract 

We provide a general formula for calculating correlators of arbitrary function of a 
Gaussian field. This work extends the standard leading-order approximation based 
on the 5N formalism to the case where truncation of the 5N at some low order does 
not yield the correct answer. As an application of this formula, we investigate 2, 
3 and 4-point functions of the primordial curvature perturbation generated in the 
massless preheating model by approximating the mapping between the curvature 
perturbation and the Gaussian field as a sum of the many spiky normal distribution 
functions as suggested by lattice calculations. We also discuss observational conse- 
quences of this case and show that trispectrum would be a key observable to search 
signature of preheating in the CMB map. It is found the forms of the curvature cor- 
relation functions for any 6N, at the leading order in the correlator of the Gaussian 
field, coincide with the standard local type ones. Within this approximation, it is 
also found that the standard formula for the non-linearity parameters given by the 
product of the derivatives of the e-folding number still holds after we replace the bare 
e-folding number appearing in the original 5N expansion with the one smoothed in 
the field space with a Gaussian window function. 



1 Introduction 



It has become a standard paradigm that the universe experienced inflation, the accelerated 
expansion, in the very early universe. During inflation, any massless scalar field acquires 
classical fluctuations on super-horizon scales originating from sub-horizon scale quantum 
fluctuations. One consequence of the inflationary scenario is that one of, or mixture of 
those scalar field perturbations convert into the primordial curvature perturbations which 
are then observed as the temperature anisotropy of the Cosmic microwave background 
(CMB) and seeds of the large scale structure (for a general review of inflation and the pri- 
mordial perturbation, for instance, see |1]). The simplest model to achieve the conversion 
is to assume that inflaton fluctuations are solely responsible for the observed curvature 
perturbations. Although such a scenario is the simplest, economical and consistent with 
the existing observational data, it is possible that more contrived scenarios as exemplified 
by the curvaton models [2H5] and the modulated reheating ones PIT] are actually realized 
in our universe. Therefore, any observable that can help us distinguish the conversion 
mechanisms deserves intensive investigation. 

In this paper, we will focus on the correlation functions of the curvature perturbation 
C, especially three-point and higher order correlation functions that describe the degree of 
non-Gaussianity of (. Higher order correlation functions have become useful to constrain 
the conversion models because of the fact that the simplest model of inflation yields only 
a negligible amount of non-Gaussianity while large non-Gaussianity up to observable level 
is the generic feature for other models (for example, see [8]). In particular, combined use 
of bispectrum and trispectrum can be powerful to constrain many models that yield large 
non-Gaussianity [9]. According to the 6N formalism [TOHTl], the curvature perturbation 
on super-horizon scales is equal to the perturbation of the local e-folding number. If 
the expansion of the universe is driven by the scalar fields, the perturbation of the e- 
folding number is produced by the perturbations of the scalar fields whose evolution can be 
calculated as if they were evolving on the homogeneous Friedmann-Lemaitre background. 
Assuming only a single scalar field sourcing ( just for simplicity, the curvature perturbation 
at point X can be written as a function of the sourcing scalar field x at the same point [15] ; 



For some models such as the curvaton models, it is a good approximation to treat the 
scalar field as the Gaussian variable. For such a case, the non-Gaussianity of ( arises due 
to the non- linear relation between x and (. In most models that fall into this category, the 
series expansion above converges so rapidly that truncation at the second order is accurate 
enough to evaluate the three-point function of (. With this approximation, we can express 
( in the form of the so-called local type as [H] 




(1) 



C — Cg + ^/NLCg) 



(2) 
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where (g = Nix is the Gaussian part of ( and /nl = is the non-hnearity parameter 
defining the strength of the non-Gaussianity of ( [15]. Current observational hmit on /nl 
by WMAP 9 year data is -3 < /nl < 77 (95% level) [H]. 

Although the use of Eq. ([2]) has been quite powerful to constrain many models that yield 
non-Gaussian curvature perturbation, there are several models for which the convergence 
of the series ([1]) is so slow that it does not admit the truncation of it at the first order for 
the calculation of the two-point function, at the second order for the calculation of the 
three-point function, etc. For instance, in a model where the curvature perturbation is 
given by trigonometric function like ( ~ cos(x/A), the expansion ([T]) becomes practically 
useless when {x^) ^ A^ (see Sec. 13.21) . Another example, which is the main subject 
of this paper, is the massless preheating model [IHllin] in which the mapping between 
X and ( can be obtained only by means of the massive numerical calculations. The 
chaotic nature of the motions of infiaton and x fi^ld which is coupled to infiaton during 
preheating and termination of the growth of field fiuctuations by the highly non-linear 
effects defeat accurate analytic derivation of the mapping. Indeed, the obtained mapping 
Cix) by numerical calculations shows that it is quite sensitive to and ( suddenly becomes 
large for some specific values of which appears as many spikes within the range of x 
given by ~ ^/{)(^■ Obviously, the expansion of the mapping given by ([1]) around some 
value of X and truncating the series at the lowest order does not work at all for obtaining 
the correct (or accurate) correlation functions of (. The latest evaluation of the mapping 
C{x) was done in [20]. In [20], the motions of the scalar fields during preheating are 
solved by performing lattice simulation, which automatically includes all the non-linear 
interactions among the scalar fields. The simulations are done for 11563 different initial 
values (i. e. values just before preheating) of x and the total number of e-fold is read for 
each simulation to obtain the mapping C(x)- Despite the resultant shape of C{x) is far 
beyond simple polynomials, the analytic fitting formula expressed by the sum of the normal 
distribution functions, each of which represents the spike, is also provided. Although they 
have discussed the possibility of generating the CMB cold spot, they have not investigated 
the higher order correlation functions of the resultant primordial curvature perturbations. 
It is then an interesting project to see what kind of forms the correlation functions of ( 
take (especially, two-point and three-point functions) as a result of the (almost) chaotic 
mapping in preheating, which motivates our paper. In Ref. [21], the authors have also 
investigated the amplitude of the primordial non-Gaussianity from preheating by using 
a smoothed quadratic function for C(x)- Om' paper differs from [21] in that we obtain 
the exact correlation functions of ( by adopting the analytic fitting formula presented 
in [20] without applying smoothing procedure from the outset. We also show that use of 
the smoothed e-folding number in the field space as the expansion coefficients of the 6N 
expansion ([1]) is correct at the leading order in the correlation function of x- 

In this paper, we first provide a basic formulation of calculating the correlation func- 
tions of Cix) [22]. By making use of this formula, we can either numerically or analytically 
evaluate the primordial non-Gaussianity in a feasible manner, for cases in which the trun- 
cation of Eq. ([T]) at the lowest order is not necessarily justified. An important assumption 
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behind our formulation is the Gaussianity of the sourcing scalar field x- After we explicitly 
verify that our formulation correctly recovers the standard local form of the three-point 
function for given by Eq. (j2]), we then apply our formulation to the models mentioned 
above such as the one where the curvature perturbation is given by a trigonometric func- 
tion and the massless preheating scenario. For this purpose, we will adopt the fitting 
formula given by [20] as the correct mapping for the massless preheating. For the physi- 
cally relevant scales, it will turn out that the three-point and the four-point functions from 
preheating still take the standard local forms. The difference from the standard one ap- 
pears in how the non-linearity parameters are related to the quantities characterizing C(x)- 
It is found that /nl from preheating is typically enormous C(IO^). We also find that the 
curvature perturbation from preheating itself does not have enough amplitude to explain 
the observed amplitude, which by necessity, requires that the total curvature perturbation 
is a mixture of dominant component and preheating component which is subdominant. 
Taking the dominant component to be the standard adiabatic Gaussian perturbation from 
inflaton, this mixture dilutes the non- Gaussianity of the total curvature perturbation and, 
as a result, /nl becomes typically (9(0.1), which is below the observational sensitivity. 
More noticeable signal appears in tnl, one of the non-linearity parameters characterizing 
the strength of a part of the four-point function [23], which is boosted typically up to 
C(IO^). These findings suggest that observational study of the trispectrum is a key to dig 
up the trace of preheating left in the curvature perturbation. 

In section |2l we provide a general formalism to evaluate the correlation function of C 
which is not necessarily written as the perturbative expansion form like Eq. (|2]). After 
demonstrating the effectiveness of the formalism by applying it to some simple models in 
section [3l detailed analysis for the case of preheating is developed in section HI In section 
m we show that, by means of the diagrammatic approach, reduction to the standard local 
type non- Gaussianity seen for the preheating case is a generic feature that holds for other 
models. The last section is conclusion. 

2 General formalism 

Let us provide a general formalism to evaluate the correlation function of the curvature 
perturbation which is not necessarily written as the perturbative expansion form [22] . Our 
primary assumption is that the curvature perturbation C, in real space is a function of a 
Gaussian scalar field x ^-t the same point: 



The second term is introduced so that {Q is zero. Without a loss of generality, we can set 
(x) = (if (x) 7^ 0; "we can redefine x by subtracting (x) from the original field.). For our 




(3) 
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#1 



purpose, it is convenient to introduce the Fourier transformation of f{x) as 

fix) = 1^ /.e^'^^ =^ f. = Jdx /(X)e-'^^ (4) 
Using fa- given above, (/(x)) becomes 

(/(X)) = J^fAewmx)a)). (5) 
For the Gaussian field x? we have the following relation, 

exp (^J £x b{x)xix)^ ^ = exp d^Xid^X2 J -^p^P^{q)e''^-^^'-'^^\xi)h{x2) 

(6) 

for arbitrary function h{x). Here, Px{k) is the power spectrum of the scalar field which 
is defined by 

{Xj:x,:.) = i27rfPxik)5{k + k'), (7) 
where Xk defined by = / e~'^'''^x{x)d^x. Using this formula, we find 

(/(X)> = / g U-"^'-, (8) 

where (x^) = {x'^{x)). Once f^ is given, we can calculate (/(x)) by performing one- 
dimensional integral. 

In a similar way, we can write down the formal expression of the A^-point function of 
/ as 

{fix,) ■ ■ ■ f{x^)) = I ^^J^^ (e^-(-)-+-+^-(-)-) 

= / (n S /..e-'^'^^^) exp (-{x') Y: -.a,e,(r.,) j , (9) 

where is defined by Cxi^ij) — {x{xi)x{x2)) / {x"^) and = \xi — Xj\. Again, once we 
know fcr, we can in principle calculate the A^-point function by performing A^- dimensional 
integral. 

Eq. ([9]) is applicable to any power spectrum of x field as long as x field is Gaussian. 
For the scale invariant power spectrum which we will assume hereafter, is given by 

p 

Px{k) = TT' -^0 : constant. (10) 

rh 



"^^Thc authors appreciate Jmi'ichi Yokoyama for suggesting this transformation. 
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Then, we have 



Here, both UV and IR cutoffs are introduced. The IR cutoff corresponds to the box size 
~ L in which perturbation is defined and UV cutoff is dependent on the model we consider. 
Carrying out the similar calculation for {x{^)x{y)) i "we finally arrive at 

-Ci{r/L) + a(g„,axr) + - ^MIh^ 

Ur) = r (12) 

In(gmax^) 

where Cj(x) is the cosine integral function defined by 

c.{x) =- r ^-^dt. (13) 



3 Application to some simple models 

Here, we present the application of our formula provided in the previous section to some 
simple models. 



3.1 Case of quadratic mapping 

For many models including curvaton model and modulated reheating model that produce 
non-Gaussian curvature perturbation, it is a good approximation to truncate the polyno- 
mial of Cix) the second order for the evaluation of the leading contribution to the three- 
point function: 

fix) = o-X + bx^, a, b : constants. (14) 

Although it is much simpler to calculate the correlation functions directly from the above 
equation as is widely done in the literature, let us here apply the formalism derived in the 
previous section just in order to verify explicitly that the formalism indeed yields the same 
results as the standard one. The Fourier transform of the quadratic expression above is 
written as 

Plugging this into Eq. ([8]) and performing the integration by parts, we find 

ifi^)) = ix') I dcr 6{a){-^aa + h)e-^^' = b{x'), (16) 
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which is obvious from direct operation of (■ ■ ■ ) onto both sides of Eq. (IT^ . In a similar 
way, the two-point function of ( becomes 

f f d (P 

(C(^)C(y)) = / daida2 S{ai)6{(r2) -ia- b- 



dai daf 

\ da2 da2 ) 

= a'{x')Ur) + 2b'{x'ni{r), (IT) 

which also coincides with the result based on the standard calculation. For the three-point 
function, after some cumbersome calculation, we arrive at 

(C(^i)C(x2)C(^3)) = 2a%{xYiUn2)Un3) + 2 perms.) + 8b' {x')%{n2)Ur23)Un3)- 

(18) 

When ( is sourced only by a single field as is considered in this paper, it is common to 
write C as a sum of the gaussian part and its square as is given by Eq. ([2]). Identifying (g 
as X aiid setting a = 1, b = I/nl in Eq. f lT8|) . we have 



(C(:?i)C(^2)C(^3)) = ^/nl((C(5^i)C(x2))(C(x2)C(5^3)) + 2 perms.) + 0{{Cf), (19) 
at the leading order in („. 



3.2 Case of trigonometric function mapping 

As the next example which is more non-trivial than the first one, let us consider a case 
where f{x) is a sine function: 

fix) = Asm {£1, (20) 

where A and A are constants. If (x^) ^ A^ is satisfied, the argument of the sine function 
becomes much smaller than unity and the Taylor expansion of the sine function up to 
the lowest order required to obtain the non- vanishing correlation function works well F^ . 
For instance, sin(x) f« x is sufficient to obtain the approximate form of the two-point 
function of (. On the other hand, for the case of (x^) ^ A^, Taylor expansion up to the 
lowest order does not provide a correct answer and it is such a case in which the formalism 
derived in this paper becomes useful. For the moment, we do not assume any magnitude 
relationship between (x^) and A^. The Fourier transform of / is given by the superposition 
of the (5-functions: 

U = A j dx sin {£ e"*^" = -niA {5{a - ^-^) - 5{a + A"^)) . (21) 

"i^^Precisely speaking, since x is 8' statistical variable obeying Gaussian distribution, it is possible that 
X takes a value much larger than A at some domain of the universe even when (x^) <C A^. But its 
probability is highly suppressed and practically it does not matter to assume x/A <^ 1 from the inception, 
that is, before taking the statistical average. 
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Using this equation, we find that {f{x)) = 0, which can be also understood from the fact 
that ( is an odd function of x- The two-point function then becomes 



#3 



(C(^)C(y)) = ^ e"^ sinh 



A2 



(22) 



We see that <C 1 is not sufficient to allow the truncation of the Taylor expansion at a 
few lowest order when (x^) ^ A^. In fact, a stronger condition, ^;^(r) <^ A^/(x^) ^ 1, 
needs to be satisfied. 

Again plugging the above /o- to Eq. and setting = 3, we can verify that the 
three-point function vanishes exactly. This should be so since the product of the three (s 
is an odd function of x- It is straightforward to extend this result to the A^-point function 
of C where A^ is any odd number and to conclude that it exactly vanishes. 

In a similar way, let us consider the case of cosine function given by 



fix) = ^cos 
The Fourier transform of / is given by 



(23) 



(24) 



Contrary to the sine function, since a cosine function is even, A^-point function of ( becomes 
non-zero even where A^ is odd number. Actually, we have 



{fix)) = Ae^ 



and 



Then, the 3-point function of ( is 

(C(^l)C(5^2)C(^3)) = 




^OS^ ' ^^x(^12) 1 cosh ( ^ex(^23; 



ix'). 



cosh 



A2 

ix') 



A2 



(25) 



(26) 



(a!) 

A2 



^""'hAru)] sinh (^4rUr23)] sinh ( ^-^^r.^] 



A2 



Un2) 



'^Hn [24], two-point function of the baryon isocurvature perturbations having the same form as Eq. ([20| 
was studied. 
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+ COsh( ^ex(^23) 



(27) 



Like a sine function mapping, the condition <^ 1 is not sufficient to perform the Taylor 
expansion up to some lowest order in this cosine mapping case, too. The above general 
expression is much complicated, but for the case with {x^)^^{r) / Is? ^ 1 we find simple 
relations given as 



(C(x)C(y)) oc Ut) 



and 



(C(X1)C(X2)C(X3)) OC ^x(^12)^x(^23)^x('^3l)- 



(29) 



This limiting case corresponds to that the curvature perturbation is a quadratic function 
of X without a linear term, that is, ( 



4 Application to massless preheating 

From the illustrations above, we have seen that the formula ([9]) can actually work to 
evaluate the correlation function of (. Now, let us focus on the case of preheating (for a 
general review of preheating, see [2S]), which is the main part of this paper. If the test 
scalar field x couples to the inflaton by the interaction: 

C = (30) 

parametric resonance generically occurs after the inflaton starts the oscillations |26f[25] . 
In particular, for inflaton potential given by the quartic so-called massless preheating 
scenario [18l[T9] , the fluctuations of the x field are not suppressed on super-horizon scales 
at the end of inflation when g^/X ~ 2 [29l[30]. Interestingly enough, the resonance band 
completely covers the super-horizon scale modes for g'^/X = 2 [19] and thus the x field 
fluctuations grow exponentially on super-horizon scales, which has invoked intensive stud- 
ies of whether the x field perturbation in the massless preheating model can generate the 
curvature perturbation on super-horizon scales during preheating [291436] (for the similar 
analysis in the other preheating models, see, for instance, [37H39] ). Linear perturbation 
analysis (linear both in metric and field perturbations) shows that the curvature perturba- 
tion ( on super-horizon scales grows exponentially at the ffist stage of preheating [291130] 



#4 



Effects of the backreaction from the field perturbations by means of the Hartree ap- 
proximation were included in [32l[33]. It was found that the backreaction terminates the 



A term "super-Hubble" instead of "super-horizon" is used in these papers because of the fact that all 
the perturbation modes we consider are actually in the causal horizon stretched significantly by inflation. 
Nevertheless, following convention widely used in the literature, we use "super-horizon" in this paper. 
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growth of (. Second order perturbation including the metric perturbations was apphed to 
the massless preheating in [IQ] and it was found that the non-Gaussianity also grows dur- 
ing preheating (although their procedure to calculate /nl was shown to be inappropriate 
in ED). 



In IH], SN formalism which takes into account the fully non-linear effects of the 
super-horizon dynamics, was applied to the massless preheating to evaluate the correspon- 
dence between ( and x (the function f{x)) for the first time. This was done by solving 
numerically the background equations of motion of the scalar fields on the Friedmann- 
Lemaitre universe. The spatial fluctuations of the scalar fields, which become important 
at the non-linear stage of preheating, are not taken into account in the analysis. It was 
observed that 6N is quite sensitive to the change of initial value of x- The origin of the 
acute sensitivity was identified as the (almost) chaotic behavior of the equations of the 
scalar fields. Thus f{x) can hardly be approximated by a quadratic expression, contrary 
to many other cases where the truncation of f{x) at quadratic or cubic order is enough to 
evaluate the bispectrum or trispectrum of the non-Gaussian curvature perturbation. The 
same conclusion was obtained in [12] in which calculations of 6N were performed for more 
number of different initial values of x- It was also suggested super-horizon scale curvature 
perturbations are generated due to the imperfect randomness of the mapping f{x)- 

Lattice calculations which automatically include the inhomogeneities of the scalar fields 
as well as the fully nonlinear interactions were implemented in [13] (see also |11]). From the 
resulting mapping /(x), /nl was estimated by fitting f{x) with the smoothed quadratic 
formula ([2]). The latest lattice calculations done by a different simulation code called 
DEFROST [l5] show that within the field range ^ (x^); ^ lot of spikes having different 
amplitudes look to appear randomly and uniformly in logarithmic interval of x- Although 
the actual form of f{x) given in [2U] is quite complex and chaotic at a first glance, the 
analytic approximation that describes the basic behavior of f{x) is proposed, which is 
given by the sum of the normal distribution: 

/(X) = E^P^^P(-^^^)' (31) 

where Xp and Hp represent the position and the width of the p-th spike, respectively. There 
are two remarks worth mentioning at this moment. First, since x is defined such that 
(x) = 0, X used in this paper is different from the one appearing in the original Lagrangian 
or the one used in [20] by a constant, which we denote by xo- This quantity represents the 
contribution from the long wavelength modes larger than the box size ~ L in which our 
perturbations are defined and effectively works at a part of the background value in our 
box. Due to the statistical properties of fluctuations, xo also varies as a stochastic variable 
as we shift the position of the box in larger space. Thus, apart from the variance (xo), 
we cannot predict a definite value of Xo and should treat Xo as a free parameter. Later, 



'^^The word SN formalism was not as popular as nowadays and a term "separate universe approach" 
was used in the paper. 
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we will study the dependence on xo of the non-Gaussianity of the curvature perturbation 
from preheating. Secondly, in the figure of f{x) given in [20], we observe two different 
types of spikes, i. e. , the spikes that have relatively large amplitudes and large widths and 
look to appear (9(10) times as we change x by one order of magnitude, and the ones that 
have much lower amplitudes and smaller widths and appear much more frequently (almost 
continuously). All the spikes in the former category have positive peaks, that is, Ap > 0, 
while the ones in the latter case have both positive and negative peaks, which appear to 
occur with approximately the equal probability. As a result, f{x) becomes positively large 
near the vicinity of any spike belonging the first category and fluctuates finely around zero 
outside them. In the following calculation, we will make an approximation that only the 
spikes in the first category contribute to Eq. ( |3T]) . Since each spike enter the correlators 
of ( in the combination ApUp, which will be shown later, the spikes in the second category 
would less contribute compared to the spikes in the first category because of the smallness 
of their amplitudes and partial offset among them. 



4.1 A mock f{x) 

Let us first consider to generate a mock f{x) that mimics the one obtained in Ref. [20]. We 
take this approach by the two reasons. The first one is that we do not have the precise raw 
numerical data of [20] but only have numbers derived from the Fig. 1 of [20] by using the 



public software l^^l which extracts the coordinates out of the graph. We will see in 14.41 that 
results using the mock f{x) show a good agreement with the ones obtained from the Fig. 1 
of [20] using the method mentioned above. The second reason, which is more fundamental, 
is that we also want to study how sensitively the slight change of /(x), i- e. amplitudes, 
positions of spikes and their widths, changes the resultant correlation functions of ( and 
to derive the generic properties of the curvature perturbation having such f{x) given by 
Eq. (13T]) . Due to many spikes whose positions apparently appear randomly in /(x), as 
was observed in [20], we decide to generate the positions of the spikes Xp randomly which 
are uniformly distributed in ln(x + xo)- Following Ref. [20], we set the range of x + Xo to 
be 0.1 = Xmin < X + Xo < Xmax = 100. Hereinafter, we present values of x in a unit of 
lO^^Afpi. The number of spikes within this interval is fixed to be 75 and the spikes are 
randomly generated within this interval. We have checked that the results are insensitive 
to both the changes of Xmin and Xmax provided that xo + \/ (x^) is well below Xmax- This 
weak dependence of the results on Xmax is due to the exponentially suppressed probability 
of realizing x + Xo > Xmax- Although this does not generically hold for the shift of Xmin, 
the weak dependence on Xmin arises due to the lower and thinner spikes for smaller x + Xo- 
We also assign the amplitude Ap by an equation Ap = a\ log]^Q(10(xp + Xo))!^'^; where 
a G (0, 4.5 X 10~^) is chosen randomly for each spike, and K,p by an equation 0.01(x + Xo)- 
These equations are employed to yield a mock /(x) that closely resembles the original one, 
apart from the fine spikes mentioned earlier. Because of the probabilistic procedure, we 
obtain different but similar /(x) by each realization. By this approach, we are able to see 



^^http:/ /www. frantz.fi/software/g3data.php 
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AN(x) 




%+Xo 



Figure 1: A mock f{x) randomly generated according to the probability distribution 
explained in the text. The parameters are such that Xo = 10 and (x^) = 10 in a unit of 
IQ-^Mpi. 



how the resulting non-linearity parameters vary by each realization and by the different 
choice of the parameters such as xo ^^'^^ ix^)- an illustration, in Fig. [H we show a 
mock /(x) for a typical case in which xo = 10 and (x^) = 10 . The result actually looks 
similar to the original f{x) provided in [2U] . 



4.2 Two-point function 

Since f{x) is given by the sum of the normal distribution, we can immediately derive the 
analytic from of which is given by 

/ — f ^1^"^ \ 
= V 27r 2^ ApKp exp ( ixp(r 1 . (32) 



Then, using Eq. ([8]), it can be easily checked that {fix)) is given by 



(33) 
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where Sp and rip are respectively the peak width and the peak position normahzed by 
a/ (x^) defined as 

ep ^ np ^ (34) 

Using the definition of C, given by Eq. ([3]), the two-point function of becomes 

{amy)) = (mm) - {f{x)f. m 

The second term on the right hand side is already given by Eq. f l33|) . The first term, by 
substituting Eq. f p2|) to Eq. Qj ior N = 2, can be written as 



Co-iCo-2 6Xp 



;(x^)(c^? + (^2) - {x{x)x{y))(^i(^2 



Pl,P2 



X exp 



^Pl ^Pl ^P2 



P2' 



1 (1 + epj< + (1 + 4)^P2 - 2ex(0^Pi^P2 



(36) 



2 (l + e^J(l + e^J-e^(r) 

Let us consider a situation where (x^) ^ Eq. ( I36l) then can be approximately written 



as 



i-e(r 



^ ^ ^pi^P2^Pi^P2 6xp 



^Pi + ^P2 ~ 2r/p,r/p,e^(r) 
2(1 -e(r)) 



(37) 



Since terms in the summation of Eq. (1371) are exponentially suppressed for |?7p| > 1, the 
summation is mostly contributed by the terms having |?7p| < 1. This means that, contrary 
to the sine function case considered in the previous subsection, -C 1 is the sufficient 
condition to perform the Taylor expansion of Eq. fl37|l in terms of and to truncate 
the series at the non-trivial lowest order that contains the non-trivial information of the 
correlator of C,. For the case of the two-point function, it turns out that we need to expand 
up to first order in ^^^^ 



where /q and /i are defined by 



#7 



-^0 — ^pCp 



e 2 



Ji = ^ ApepTipC 



(3^ 



(39) 



Applying again the approximation (x^) ^ to Eq. fl5^ . we find 

(/(X)) = h. 



(40) 



Although not written explicitly, both /q and /i depend on xo- 
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Figure 2: Two-point function of C{x) for realization (Jq = 0.67 x 10~^, Ji = 0.19 x 
10~^, I2 = 0.46 X 10"^). The solid line represents the exact correlation function calculated 
from Eqs. (!33l) and (136|) . The dashed line represents the approximated correlation function 
given by Eq. (14T1) . In both cases, the correlation functions are decreasing function. At the 
left edge where r/L = 10~^°, we find ~ 0.38 and it decreases to ^-^ ~ 0.017 at r/L = 1. 

It is clear that the first constant term in Eq. fl38l) exactly cancels with the second term in 
Eq. fl35|) . In the diagrammatic language, the first term and the second term in Eq. f l38|) 
correspond to the disconnected and the connected diagram, respectively, and the non- 
trivial information are contained in the connected diagram. Then two-point function of ( 
becomes 

{ax)cm = iiUr) + o{e^). (4i) 

This shows that in the regime where .^^^ ^ 1, the two-point function of ( is proportional 
to and the scale dependence of C,^ is stored in {({x)({y)). 

Fig. [2] shows the two-point function of C{x) for the same parameter choice, Xo = 
10) (x^) = 10. The solid line is the exact correlation function calculated from Eqs. f l55|) 
and (1361) . The dashed line represents the approximated correlation function given by 
Eq. ( 14T|) in which second and higher order terms in are dropped. Here, we have assumed 
Qma^L = 10^^ which appears in ^xi''^)- "^^^^ is because gmax is roughly equal to the Hubble 
scale at the end of inflation for the massless preheating case with g"^/ A = 2 and taking L 
to be the size of the present Universe we have QmaxL ~ 10^^ (to^gcv)' At r/L = 10~^°, 

becomes as large as 0.38 and it monotonically decreases to 0.017 at r/L = 1. The 
current observations can probe the primordial fluctuation from the current Horizon scale 
down to the several order of magnitude smaller than it. This leads to the huge hierarchy 
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L >r ^ Q'max- For such an observationally relevant range, C.xi'f') given by Eq. (fT2|) becomes 
mildly smaller than unity (typically 0.05 ~ 0.1). Although is suppressed only mildly in 
this case, in order to capture the basic point of our result in analytic way, we regard as 
the expansion parameter and Taylor-expand the correlators with respect to it as can be 
seen in Fig. El 

It is interesting to consider the meaning of the first term on the R. H. S. of Eq. fHT]) . 
Using the definition of Ji given by Eq. fl39|) . we have 



^1 = IZ^p'^pXpe"^ = -7=-^ / dx' fix + Xo)x'e~^. (42) 
\X } p V27r(x2) J 

We now introduce a Gaussian window function Wg{x) defined by 

W^.(x) = ^=i==e-A. (43) 

This satisfies the normalization condition J dx Wg{x) = 1 and effectively cuts off con- 
tributions from \x\ ^ a/ (x^)- With the Gaussian window function, Ii can be written 

as 

h = / dx' fix' + Xo)x'Wg{x'). (44) 

On the other hand, we can define f^{x) @ which is defined as a smoothed f{x) iii the 
region (x - V (x^)' X + a/(x^)) by 

/''(X) = I dx' f{x + x')Wg{x'). (45) 
We can verify that Ji is related to the first derivative of f^{x) evaluated at x = Xo; namely, 

(46) 



dx 

Therefore, we can rewrite Eq. (14T]) as 



x=xo 



2 



{mem ^ y^^j (x(^)x(y)). (47) 

This result shows that even when the original e-folding number is spiky and cannot be 
approximated as linear expression of the scalar field at all, at the leading order in the 
two-point function of C, is given by the first derivative of the e-folding number, just like 
the standard 5N formula, provided it is smoothed over the range ~ a/ (x^)- As we will 
see later, this correspondence holds for the three-point function and even for the A^-point 
function of a general function /(x) (at leading order in ^^). 



This quantity was computed numerically in f20[|42| to estimate the resulting curvature perturbation 
on super-horizon scales. 
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4.3 Non-Gaussianity 

It is a straightforward and standard calculation to perform the three-dimensional integral 
of Eq. dH]) for f^r given by Eq. ( 1321) . However, since the resulting expression is complex and 
long, we write down the expression after applying the approximation (x^) ^ as what 
we did for the case of the two-point function. The result is given by 

Pl,P2,P3 

where S and T are defined by 

^ = 1 - ^^(^2) - e?(r23) - e^(r3i) + 2e^(ri2)ex(^23)ex(^3i), (49) 
^ = ^1 + ^2 + ^3 - 2{r7ir72(^x('^i2) - '^x(^3i)^x('^23)) + 2 perms.} 

-vl^Ur,,) - rie^^) - Vie^in,), (50) 

where = |xj — In the regime where <^ 1, the three-point function of ( is 

(C(xi)C(^2)C(^3)) = {h - h)ll (ex(^i2)ex(n3) + 2 perms.) + 0(g). 

^ ^^72^ ((C(xi)C(5^2))(C(5^i)C(x3)) + 2 perms.) , (51) 
-'i 

where I2 is defined by 

/2 = $^Apep^Je-^. (52) 

p 

Interestingly, the last equation takes the same form as Eq. f|T9|) . i. e. , the three-point 
function of C, at physically relevant scales becomes the standard local type one common to 
many models such as curvaton model and the modulated reheating one. The corresponding 
/nl is given by 

-/NL = ^2 • (53) 

In the same way as we derived Eq. (jH]), we can write Jq and I2 in terms of the Gaussian 
window function as 

dx fix + Xo)W,{x'), ^2 = ^1 dx' fix + Xo)x"W,ix'). (54) 

We can show that the second derivative of f^ix) defined by Eq. (H5l) is given by 

d'f^iXo) 1 



dx^ (X^) 
Then, we have 

6 . fxx 



ih-h). (55) 



5/- - (56) 
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Therefore, /nl is also given by the standard expression based on the SN formah sm prov ided 
the e-folding number is understood as the smoothed one over the region ~ a/7x^0- 

Given the current situation where cosmic measurements turn out to be useful to probe 
the four-point function, or trispectrum (for instance, see [IS]), it is intriguing to go one 
step further from the three point function and to evaluate the four-point function. The 
calculation is tedious but straightforward. Again with the approximations that (x^) ^ i^p 
and .^1^, ^ 1, we find that the four-point function reduces to 

(C(^i)C(x2)C(^3)C(^4)) ^ 7/°^ ((C(^i)C(^2))(C(^2)C(^4))(C(^3)C(^4)) + H perms.) 



ih - 3/ 



If 



- ((C(xi)C(x2))(C(^i)C(x3))(C(^i)C(^4)) + 3 pern<Sr) 



where 



p^pvle 2 . (58) 



We see that the four-point function is contributed by two distinct parts. Indeed, this is ex- 
actly the same form as the standard local type four-point function which is parametrically 
represented as 

(C(^l)C(^2)C(^3)C(^4)) = T-NL ((C(^l)C(^2))(C(^2)C(^4))(C(^3)C(^4)) + H permS.) 

54 

((C(^i)C(x2))(C(^i)C(^3))(C(^i)C(^4)) + 3 perms.|59) 

where tnl and (^nl are non-linearity parameters for the four-point function [23]. Comparing 
this with Eq. ( 157|) . we can read 

_{h-Iof_36 , 54 _ (/3-3/1) 
Tnl ji ^/nl' 25^^^^ ~ If ' 

Our results show that the non-linearity parameters are given by the combination of the 
moments = 0, 1, 2, 3). Each moment can be numerically evaluated provided position rjp, 
amplitude Ap and width ep of each spike are known, which we will carry out in a following 
way. But before going to the numerical results, we can make a crude numerical estimate of 
the non-linearity parameters. From the figure in [2U], for the relevant range of x + Xo? we 
find that Ap = 0(10"^) - 0(10"^) and ep = O{10-^) - O(O.l) and A^, = C(10) spikes. As 
is clear from the expression of Jq (see Eq. (l39l) ). only the spikes having |?7p| < 1 contribute 
to Jo and we have Jq ~ {Ap){ep)Ns ~ 0(10"^). Similarly, we expect h ~ 0(10-^). On 
the other hand, since terms in the summation of Ji are either positive or negative with 
the (almost) equal probability, they may partially cancel out each other. As a result, 
we would expect |Ji| is at most O{10~^). Therefore, unless accidental cancellation to 



By adopting values of and provided by [50], this fact was used in [5T] to compute /nl from 
the massless preheating. 
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Figure 3: Distribution of the non-linearity parameters /nl and (^nl for the case Xo = 10 
and (x^) = 10 in a unit of 10~^Mpi. The number of reahzations is 200. Since both /nl 
and (7nl can be both positive and negative and scatter over a several orders of magnitude, 
both axes are presented logarithmically. Each quadrant is defined for (/nliS'nl) plane. 
For instance, the 2nd quadrant corresponds to a region in which /nl < 0, (/nl > 0. 

have I/2 — ^ 1 occurs, as a crude estimation, we have |/nl| ^ O(IO^). Similar 

consideration yields \g^i\ ^ 0{1Q^^). As we will see below, this estimation indeed provides 
the correct orders of magnitude of the non-linearity parameters. 

With the same parameters, we performed 200 realizations of the mock /(x) and calcu- 
lated /nl and (/nl for each realization. Fig. [3] shows the resultant distribution of the non- 
linearity parameters plotted in logarithmic scale for fixed parameters xo = 10, (x^) = 10. 
Since both /nl and (/nl can be positive or negative, the distributions are shown for each 
quadrant defined for (/nl^S'nl) plane. We see that the distribution is highly concentrated 
in the region specified by 10"^ x /nl = C(0.1)-C(10) and lO"^^ ^ ^^^^ = C(O.l) -C(IOO), 
which roughly agrees with our earlier analytic estimation. Also, we see a clear tendency 
that larger |/nl| entails larger Is'nlI- We also find that the distribution is wide and covers 
a several orders of magnitude. The large |/nl| is realized when Ii becomes much tinier 
compared to Jq and I2 due to the accidental cancellation among terms in the summation 
given by Eq. fl39|l . Since both /nl and (/nl have Ii in their denominators, (/nl generically 
becomes large as well when large |/nl| is realized. Actually, the distribution for higher 
|/nl| lies on a line |(/nl| = I/nlI^^^, which supports the explanation due to the cancellation 
of Ii. The opposite case in which IJ2 — /o|//o ^ 1 accidentally happens is also possible and 
/nl becomes small in such a case. Since the numerator of (/nlj which is given by J3 — 3Ji, 
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Figure 4: The function f{x) for one realiza- Figure 5: The function f{x) for one reahza- 
tion that yields lO"*' x /nl = 0.0055. tion that yields 10"^ x /nl = 1.1 x 10^ 

is independent of that of /nl, the vanishing of IJ2 — /o| does not result in the vanishing 
of I3 — 3/1 in general, which explains the disappearance of the positive correlation that 
exists for higher |/nl|- This consideration suggests that, in principle, either of any very 
large |/nl| (and Is'nlI) or very small /nl can be realized after we perform sufficiently large 
number of realizations. In other words, we can say that the non-linearity parameters are 
quite sensitive to the realizations of /(x)- This would imply that the (slightly) imperfect 
modelling of f{x) has a possibility of leading to a significant error on the estimation of the 
resulting non-linearity parameters. Just for illustration, in Figs. Il]and[5l we show two dif- 
ferent fix) due to the different realizations with the same parameters xo = 10, (x^) = 10. 
These correspond to the extreme cases where 10~^/nl = 0.0055 for the former case and 
10~^/nl = 1-1 X 10^ for the latter. It is easy to recognise that they have the similar feature 
but different in detail, but not easy, without implementing the numerical calculations, to 
say if they produce quite different /nl or not. 

Very large /nl corresponding to right-hand region of Fig. 12] should be taken with caveat. 
As we mentioned earlier, such very huge /nl is realized when /i becomes tiny. In such 
a case, the higher order terms in appearing in the correlation functions, which are 
neglected as being subdominant, are no longer suppressed. In order to clarify this point 
more quantitatively, let us extend the expansion of Eqs. (14T|) and (I5T]) to one more higher 
order; 

(amy)) = iiur) + \{h - hfe^{r) + o{ii), m 

(C(^i)C(^2)C(^3)) = {h - {Un2)W>^vi) + 2 perms.) 

+hi{h - 3h){h - /o)(ex(n2)e'(r23) + 5 perms.) 

+ ih - IofUn2)Ur23)Uri3) + 0(g). (62) 



18 



We find from tlie first equation that the higher order terms become important if we 
have If < |(/2 — -^o)^^x — -'-^"^^^x- terms of /nl, we can rewrite this condition as 
/nl ^ Therefore, if /nl becomes larger than 10^/^^, the standard expression of 

the three-point function written as the two product of the two-point functions, which is 
given by Eq. (fT9l) . loses its vahdity and the use of the constant /nl to characterize the 
strength of non-Gaussianity becomes pointless. In stead of the standard relation, in the 
limit Ji — 7- 0, we have 

(C(^l)C(^2)C(^3))' = 8(C(fl)C(^2))(C(^2)C(^3))(C(^3)C(^l)). (63) 

Thus the three-point function is completely determined by the two-point function without 
unknown parameter. 

We also study how the non-linearity parameters vary as we shift xo fo^^ ^ fixed f{x) 
generated randomly (one realization). Fig. [6] shows the result in which xo is varied from 
1 to 50 with the condition (x^) = 10. As we can see, /nl becomes quite tiny (and even 
becomes zero) at some specific values of Xo- 9nl also does, but at different values of Xo? 
which is expected because of the reasoning we made earlier. We also observe that /nl is 
significantly enhanced (and even diverge) at some specific values of Xo for which (^nl is also 
amplified as well. This can be again understood as the accidental vanishing of Ii due to the 
cancellation of terms entering the summation of Ji. Although this is for one realization, 
we checked that the qualitative feature remains the same for other realizations. Thus, it 
is a general consequence that we have /nl = C(IO^) and (7nl = C(10^^) for most of the 
values of xo, but they happen to become either tiny or huge if xo is suitably fine-tuned. 

To summarize our findings, we have two sources of the uncertainties that deaden the 
predictability of the model on the non-linearity parameters. The first one is due to the 
incomplete imitation of the mock /(x)- As we have seen, two mock /(x) that look similar 
occasionally yield quite different values of /nl and (^nl with some small probability. This 
issue may be, in principle, resolved by knowing the correct form of the real /(x) precisely, 
which is beyond the scope of this paper. Furthermore, even if we know the correct form 
of f{x), it is not clear of whether the Gaussian fitting given by Eq. ([SI]), which is essential 
in the sense that it allows analytic study of the correlation functions to some extent, is 
approximate enough to derive the reliable values of the non-linearity parameters. This is 
a difficult problem and we do not pursue it in this paper. Instead, we take a conservative 
stance that our approach can only estimate the likely values of /nl and (^nl- 

The second uncertainty is due to the indeterminacy of Xo, which is more fundamental 
than the first one. Since Xo is a statistical quantity reflecting the large wavelength modes 
outside the box size, it is impossible in principle to predict a definite value of Xo- As a 
result, we cannot predict the values of /nl and (^nl, although it may be allowed to say at 
least that /nl = C(IO^) and (^nl = C(10^^) are likely to happen. 
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Figure 6: Dependence of the non-linearity parameters on xo for one realization with 
(x^) = 10 in a unit of 10~'^Mpi. The vertical axis represents 10"^ |/nl| and 10~^^ IQ'nlI in 
logarithmic scale. 

4.4 Consistency check of the mock mapping 

So far, our study is done with the mock mappings which are not the ones obtained by 
the lattice calculations. It is therefore interesting to compare the results based on the 
mock mappings with the ones based on the lattice calculations. Since we do not have the 
original numerical data of f{x), we have read off the coordinates of f{x) out of the Fig. 1 
presented in [20] by using the public software. In order to extract the position, height and 
width of each spike, we smoothed the obtained f{x) by the Gaussian window function 
with a width 0.015(xo + x) which is smaller than \/J^x^ but is large enough to eliminate 
the fine spikes. Fig. [7] shows the resultant smoothed graph. Height of each spike is a little 
bit smaller than the original one given in [20], which is caused by the smoothing. The 
width is, instead, more broadened to compensate the loss of the height. From this data, 
we extract the position, height and width of each spike. We set the minimum height of 
spike we count in the analysis to be lO^AA^. We choose this value to avoid the inclusion 
of small spikes which appear to deviate from normal distribution function. 

Fixing Xo and {x'^), we can calculate the moments Iq, /i, ■ ■ ■ by using the Eqs. fl39|) and 
so on. Fig. [H] shows the dependence of /nl on xo with (x^) = 10 by using the formula (jHSD- 
We see that Fig. |H] looks quite similar to Fig. IH] in a sense that /nl is typically C(IO^) but 
accidentally becomes very small or large for some particular values of xo- This supports 
that our findings derived in the previous section by using the mock mappings capture the 
essential part and are not modified significantly by the use of the real mapping. 
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Figure 8: Dependence of /nl on xo with 
Figure 7: Smoothed /(x) of Fig. 1 presented {x'^) = 10 in a unit of 10"'^Mpi. The ver- 
in [20]. Unit of the horizontal axis is IQ-^Mpi. tical axis represents 10"*^ |/nl| in logarithmic 

scale. 



4.5 Observational consequences 

So far, our analysis is based on the assumption that the total curvature perturbation orig- 
inates from preheating. This cannot describe the primordial perturbation of our Universe 
since the produced curvature perturbation from preheating does not have enough ampli- 
tude to explain the observed amplitude. Assuming Ji = 10"^ and = 0.04 corresponding 
to r/L = 10"^ as typical values, we have {C{x)({y}) = 4 x 10"-^^ from Eq. fHTl) . At the 
Sachs- Wolfe region corresponding to low i in the CMB multipole i < 60, the curvature per- 
turbation is related to the CMB temperature anisotropy by ({x) = — 5AT(x)/Tcmb, where 
^CMB ~ 2.7 K is the CMB mean temperature. Using the WMAP observations which detect 
^ 10^ /iK^ of the power of the temperature fluctuations [IT], we find {C{x)({y)) ~ 7 x 10 
is required from observations. This is about four orders of magnitude larger than that from 
preheating, which is observed in Fig. [2l Therefore, we need additional perturbations that 
account for the dominant part of the total curvature perturbations. Let us assume that 
the standard adiabatic perturbation from inflation plays this role. With this assumption, 
the total curvature perturbation is a mixture of two contributions, which can be written 
as 

C = Cinf + Cpre. (64) 

It should be now understood that all the calculations done in the previous subsections 
are for (pj-^. It is known that the single field inflation with the canonical kinetic term 
predicts deviation of Cinf from Gaussianity is suppressed by the slow-roll parameters. If 
the Lagrangian of the inflaton field is more complex or there are multiple fields responsible 
for inflation, non-Gaussianity of Cinf can be boosted up to observable level in general (see 
for instance [8] and references therein.). In this section, we assume that Cinf is Gaussian 
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just for simplicity. From the requirement that ( explains the observed amplitude of the 
CMB temperature anisotropy, we have 



(amy)) ^ (Cinf(x)Ci„f(y)). (65) 

Therefore, to fit this with the observed slope of the two-point function (spectral index) 
is achieved simply by picking up the suitable inflaton potential. Dynamics of preheating 
does not enter this game. On the other hand, from the Gaussianity ansatz for d^f, we also 
have 

(C(^l)C(^2)C(^3)) ~ (Cprc(^l)Cprc(^2)Cprc(^3))- (66) 

In the region where the lowest order approximation given by Eq. ( 15 ip works well, Eq. (16 6 p 
can be written as 

(C(^i)C(^2)C(^3)) ^ s'^-^^ ((C(xi)C(x2))(C(i^i)C(i^3)) + 2 perms.) , (67) 



where s is defined by 



^ (Cpre(x)Cprc(y)) , ^ 

(amy)) ' ^ ^ 

and represents the relative contribution of the curvature perturbation of the preheating 
origin to the total two-point function. By definition, we have s < 1. Strictly speaking, s 
is a function of \x — y\ (although its dependence on r is logarithmic) and the factorization 
of Eq. ( l68l) is mathematically inconsistent. Nonetheless, this factorization is still useful 
to estimate the magnitude of the three-point function if it is observationally important or 
not. Taking (CprcCprc) ~ 4 x 10""^^ and {(() ^ 7 x 10~^^ as representative values, a typical 
value of s is estimated as s = 6 x 10~^. Comparison of Eq. ( 167|) with Eq. ( !T9|) enables us 
to define the effective /nl given by 

^/nl = s ^2 • (69) 

Choosing s to be 6 x 10~^, we see that the effective /nl gets smaller than the original 
one, for instance, shown in Figs. Eland El by a small number = 4 x 10~^. Interestingly 
enough, while the original /nl is too large to be compatible with the existing observational 
constraint |/nl| ^ 100 except for the cases in which xo is fine-tuned to some specific 
values, the suppression by the factor makes the effective /nl to be typically in the range 



(9(10~^) — (9(0.1) l^lj, which is below the detectable amplitude by observations. One may 
wonder if one can boost /nl to much higher value than 0(0.1) by making Ji be accidentally 
small. However, this is not so since /nl in the mixed case is proportional to rather than 



'^^^ The general relativistic second order effects produce /nl = Such effects are not considered in 

this paper. 
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to If ^. Thus, smaller Ji results in smaller /nl and /nl ^ C'(l) is the robust upper bound 
on /nl in the massless preheating of our case. 

In a similar way, the four-point function can be written as 

(C(fl)C(^2)C(^3)C(^4)) ~ (Cprc(Xl)Cpre(x2)Cpre(^3)Cpre(x4)) 



n 



4 



s^((C(xi)C(x2))(C(x2)C(x4))(C(^3)C(5^4)) + 11 perms.) 



^(4_3^^3 ((C(fi)C(x2))(C(^i)C(x3))(C(i^i)C(^4)) + 3 pern<§:^) 



Correspondingly, the non-linearity parameters are given by 



ih - hf „3 36 2 54 (/s - 3/i) ^3 

TNL = = ^/nL, ^^7NL = • (71) 



Thus, tnl is enhanced by s = (9(10 ) compared with /^l- Actually, tnl is always 
larger than or equal to irrespective of the underlying model [17H19]. Particularly, 

as the above equation shows, enhancement of tnl compared to /^l is the general feature 
for the case where the curvature perturbation is a mixture of the dominating Gaussian 
perturbation and the sub-dominant non-Gaussian perturbation [nOlEI]- In light of the 
situation that Planck satellite can constrain tnl to be up to ~ 600 [52], observational 
studies of trispectrum can be quite useful to examine if the trace of preheating is left in 
the non-Gaussianity of the curvature perturbation. The second equation of ( ffTl) shows 
that we may expect (^nl to be 0{1). This is a several orders of magnitude smaller than 
the sensitivity expected to be accomplished in future observations [53], but is still larger 
by a few orders of magnitude than the one coming from the standard adiabatic fluctuation 

Cinf- 



5 General cases 

Here, we have seen that, for the curvature perturbations sourced by preheating, three- 
point and four-point functions at the leading order in C,x take exactly the same forms as 
those of the standard local type non-Gaussianity for which the curvature perturbation is 
given by the expansion of x as 

C(^) = 5" ^iVnX"(^) - average, (72) 

n=0 

at the leading order in C,x- In this section, we will show that this is not the special 
property limited to the curvature perturbations from preheating, but the generic feature 
of any curvature perturbation given by Eq. ([3]), for which expansion in terms of x needs 
not to be a good approximation, provided x obeys Gaussian statistics. Thus, the only 
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assumption we will make in this section is that is small enough to allow us to pick up 
the leading terms in any correlator of (. We will thus leave the function form of f{x) 
unspecified. 

In the present situation, we can Taylor-expand Eq. iQ in terms of C,x truncate it at the 
leading order. Thus the resultant expression consists of terms, each of which contains the 
product of Cxi^ij)- Each fully connected term in the correlator, which we are interested 
in, corresponds to the product such that any (z = 1, ■ ■ ■ , A^) appears at least once in 
the argument of and a corresponding diagram, which has A^-vertices to each of which 
Xi is attached without duplication and is constructed by drawing line between Xi and xj 
if we have Vij in is simply connected. For instance, ^^(^is)^^^'^^^) appearing in the 
four-point function does not give the connected term, while ^x('^i2)^x('^23)^x('^34) does. At 
the leading order in corresponding connected diagrams are tree diagrams, which do not 
have any loop inside, because any connected loop diagram is made by attaching lines to 
some connected tree diagram with the same number of vertices, which is accompanied by 
additional powers of It turns out that any connected tree diagram having A^- vertices 
has A^ — 1 internal lines. This means that any leading connected term in the A^-point 
function of ( contains a product of (A^ — 1) C,x 's- Therefore, at the leading order in we 
have 



{as^)■■■ax^)) = {-{x')f~'c 



N-1 



n ^xirp.,. 



p=i 




N 



nd(Ti _ 
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(73) 



where C[- ■ ■] means we take only connected parts of [■ ■ ■]. As mentioned above, we can 
associate each term in Eq. f l75]) with the corresponding connected diagram. Since all 
the possible combinations (■ ■ ■ ,Pi,Pj, ■ ■ ■ ) in the product of 's appear in Eq. f l75]) . all 
the possible connected tree diagrams contribute to Eq. (1731) . Since the connected tree 
diagrams having A^-vertices that are not isomorphic to each other show different scale 
dependence due to different appearance of the combination of the product of C,x{^Pi,Pj)- 
The number of connected tree diagrams is equal to that of the independent parameters 
required to parametrize the A^-point function. This is exactly the same property for the 
A^-point function of the standard local type curvature perturbation of the form f l72l) at the 
tree level [MIES]. Therefore, Eq. fl73|l has exactly the same structure as the one for the 
standard local type curvature perturbation. 

Conversely, from a given connected tree diagram, we can construct a corresponding 
non-linearity parameter contributing to Eq. (1751) according to a following rule. For a vertex 

to which a vector Xi is attached, we assign — %;^/o-ie~"^°"»^crf% where gi is a number of lines 
coming out of the vertex. Without a loss of generality, we can assume gi < g2 ■ ■ ■ ^ Qn- 
For a line connecting the vertices Xi and Xj, we assign ^xi'^^ij)- then multiply all of 
them and integrate it over (cxi, ■ ■ ■ ,crAr), which gives one term in Eq. f l73p . In order to 



make things simpler, let us define the moments J„ 



n 



0,1, 



by 



'J n 



da 
2^ 



fae~ 



a 



(74) 
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A connected tree diagram having A^- vertices can be characterized by a sequence consisting 
of integers ((?!, ■ ■ ■ ,gN)- Therefore, a non-hnearity parameter f{gi,--- ,gN) associated 
with the corresponding diagram can be written as 



TV 



f{9u--- ,9N) = J^''^''-'''l['^,.. (75) 

1=1 

As illustrations, let us consider the three-point function. There is only one connected tree 
diagram with three vertices characterized by {gi, g2, gs) = (1, 1, 2). Then, the non-linearity 
parameter can be written as 

^/nl = JfVs. (76) 
5 

It can be explicitly verified that this equation exactly recovers Eq. fl53l) for the preheating 
case. In a similar way, for the four-point function, there are two connected tree diagrams 
characterized by (1, 1, 2, 2) and (1,1,1,3). The former case gives tnl given by 

36 

Tnl = Ji^Jl = ^/nl- (77) 



The latter case give (7nl given by 



54 

—g^L = Ji^Jz- (7^ 



This also coincides with Eq. f l60p for the preheating case. 

Finally, let us show that just as in the case for massless preheating in the last section, J„ 
is related to the n-th derivative of the smoothed e-folding number. The Fourier transform 
of the Gaussian window function defined by Eq. ( 143|) can be written as 

Wg{a) = jdx W,{x)e-'^'' = e'^'^'. (79) 

With this Fourier-transformed window function, we can write the smoothed e-folding 
number f^{x) as 



Thus, its n-th derivative becomes 



■n rR 



where we have used Eq. ( TMl) in the last equation. This shows that, to compute the 
correlation function of C, we can apply the standard diagrammatic rule following from the 
5N formalism, namely to assign ra-th derivative of the e-folding number to the vertex where 
n lines are attached, simply by replacing the bare e-folding number with the smoothed 
e-folding number. 
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6 Conclusion 



There are some inflationary models where the resultant curvature perturbation ( is a. func- 
tion of the Gaussian scalar field but its dependence cannot be approximated by quadratic 
expression. Representative models include the trigonometric mapping and the massless 
preheating models. In this paper, we provided a general formulation to calculate the 
correlation functions of the curvature perturbation which is an arbitrary function of the 
Gaussian scalar field x- The iV-point function of the curvature perturbation is written as 
A^-dimensional integral, which can be evaluated once the Fourier transform of the mapping 
between ( and x is given. We applied the formalism first to the common local type 
non-Gaussian curvature perturbation and verified that it reproduces the standard local 
form of the three-point function. We then considered more non-trivial cases including the 
model in which the curvature perturbation is the trigonometric function of the scalar field. 
Due to the appearance of the Dirac's 5-function, this model allows analytic evaluation of 
the correlation function. 

By using the fitting formula provided in [20] which approximates each spike of the 
mapping C{x) as the normal distribution function, we could analytically perform the inte- 
grals for the two-point, three-point and four-point functions of (. For physically relevant 
scales, forms of the correlation functions reduce to the standard local type ones. In par- 
ticular, while the non-linearity parameters defining the strength of the non-Gaussianity of 
( are related to the derivatives of the e-folding number with respect to x the standard 
case, the role of the derivatives is replaced by the moments defined as sums over spikes 
with suitable weights for the case of preheating, or equivalently, by the derivatives of the 
e-folding number smoothed in the field space. Due to the lack of the knowledge of the 
precise positions of the spikes and their widths and amplitudes, we randomly generated 
the mapping of mock C{x) 's having the similar properties to the original one found in [20] 
and studied how the resultant non-linearity parameters differ according to different realiza- 
tions. We found that the non-linearity parameters take naively estimated values for most 
of the realizations but some realizations yield very large or small values due to occasional 
cancellation among contributions from each spike. Even after the parameters appearing 
in the original Lagrangian are fixed, statistical indeterminacy of the average value of x iii 
our observable universe also causes another uncertainty of the non-linearity parameters. 
We therefore studied dependence of the non-linearity parameters on the average value of 
X for a fixed realization. It was found that for some specific values (these values differ 
by different realizations) the non-linearity parameters are sharply enhanced or diminished 
compared to the naively estimated values. These analyses show that we can only estimate 
the likely values (order of magnitude) of the non-linearity parameters. 

Due to the insufficient amplitude of the curvature perturbation from preheating, the 
curvature perturbation in our observable universe, if the massless preheating indeed hap- 
pened in the early universe, must be a mixture of another dominant component and the 
subdominant one of preheating origin. Assuming the dominant component to be the stan- 
dard adiabatic Gaussian perturbation coming from infiaton, the mixture reduces the non- 
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Gaussianity of the total curvature perturbation to the observationally allowed range. The 
trispectrum is relatively amplified compared to the bispectrum and use of both observables 
will be useful to search preheating signature buried in the curvature perturbation. 

By a diagrammatic approach, we also showed that when the correlation function of 
X can be treated perturbatively, the forms of the correlation functions of C coincide with 
the standard local type ones at the leading order approximation. While the non-linearity 
parameters are given by the product of the derivatives of the e-folding number in the 
standard case, they are given exactly the same manner even in our case provided we re- 
place the bare e-folding number with the one smoothed in the field space with a Gaussian 
window function. Diagrammatic rule for converting the diagram to the corresponding 
non-linearity parameter was also provided. 
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